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There have been extensive experimental search for possible exciton superfluid in semiconductor 
electron-hole bilayer systems below liquid Helium temperature. However, exciton superfluid are 
meta-stable and will eventually decay through emitting photons. Here we show that the light emitted 
from the excitonic superfluid has unique and unusual features not shared by any other atomic or 
condensed matter systems. We evaluate angle resolved photon spectrum, momentum distribution 
curve, energy distribution curve and quasiparticle excitation spectrum in the exciton superfluid and 
comment on relevant experimental data in both exciton and exciton-polariton systems. 



1. Introduction. An exciton is a bound state of 
an electron and a hole. Exciton condensate was first pro- 
posed more than 3 decades ago as a possible ordered state 
in solids 0,0, 0|- But so far, no exciton superfluid phase 
has been observed in any bulk semi-conductors. Recently, 
degenerate exciton systems have been produced by differ- 
ent experimental groups with different methods in quasi- 
two-dimensional semiconductor GaAs/AlGaAs coupled 
quantum wells structure 0, S S When the dis- 
tance between the two quantum wells is sufficiently small, 
an electron in the conduction band in one quantum well 
and a hole in the valence band in the other quantum well 
could pair to form an indirect exciton which behaves as a 
boson in dilute exciton density limit. It was established 
that the indirect excitons in EHBL has at least the fol- 
lowing advantages over the excitons in the bulk: (1) Due 
to the space separation of electrons and holes, the life- 
time T ex of the indirect excitons is 10 3 ~ 10 5 longer than 
that of direct ones, now it can be made as long as mi- 
croseconds. (2) Due to the relaxation of the momentum 
conservation along the z direction, the thermal lattice 
relaxation time tt, of the indirect excitons can be made 
as 10 ~ 3 that of bulk excitons, so r ex ^> T£ is well satis- 
fied. (3) The repulsive dipole-dipole interaction is crucial 
to stabilize the excitonic superfluid against the compet- 
ing phases such as bi-exciton formation and electron-hole 
plasma phase. So EHBL is a very promising system to 
observe BEC of in-direct excitons. The quantum degen- 
eracy temperature of a two dimensional excitonic super- 
fluid (ESF) can be estimated to be T% x ~ 3K for ex- 
citon density n ~ 10 10 cto -2 and effective exciton mass 
to ~ 0.22mo where too is the bare mass of an electron, 
so it can be reached easily by He refrigerator. Indeed, 
as temperature is decreased from ~ 2QK to ~ 1.7 K, the 
spatially and spectrally resolved PL peak density center- 
ing around the gap E g ~ 1.545eV increases the ex- 
citon cloud size decreases to L ~ 30/xm, the peak width 
shrinking to ~ lmeV at the lowest temperature ~ 1.7 K. 
All these facts indicate a possible formation of exciton 



condensate around 1.7 K . 

In this paper, we will study quantum nature of photons 
emitted from the excitonic superfluid phase in semicon- 
ductor electron-hole bilayer systems. We comment on 
current PL experimental data and also propose possible 
future experiments such as angle resolved power spec- 
trum to test the existence of exciton condensate in the 
EHBL system. The phase sensitive homodyne measure- 
ment to detect the two mode squeezing spectrum and 
HanburyBrown-Twiss type of experiments to detect two 
photon correlations and photon statistics will be pre- 
sented in a separate publication I* 1 this paper, for 
simplicity, we ignore the spins of excitons, the effects of a 
trap and disorders. Their effects are important and will 
be investigated in separate publications. 

2. Exciton-Photon Hamiltonian. The total 
Hamiltonian in grand canonical ensemble is the sum of 
excitonic superfluid part, photon part and the coupling 
between the two parts H t — H — fxNt — H s t + H p h + Hi nt 
where : 

H sf = ^{Et-^b\b i+ ^VMb\J p+ ^ k 

k kpq 
k 

H int = ^2[ig(k)a k bl + h.c.}. (1) 

k 

where A is the area of the EHBL, the exciton en- 
ergy Ej. x = k 2 /2M + Eg — Eb, the photon frequency 

k>fc = v g\Jk 2 + k 2 where v g = c/^/e with c the light 
speed in the vacuum and e ~ 12 the dielectric constant 
of GaAl, k = (k, k z ) is the 3 dimensional momentum, 
Vd(q) = ^f^(l — e~ qd ) is the dipole-dipole interaction 
between the excitons [T^, V^flf] ^> d) = e 2 d 2 / \r\ 3 and 
Vd(q = 0) = 27Te e d leads to a capacitive term for the 

density fluctuation The g(k) ~ tk\ • Dk x L^ 1 ^ 2 is 
the coupling between the exciton and the photons where 



2 






o k 
(c) 



FIG. 1: The energy spectrum and the decay rate of the exciton 
versus in-plane momentum k. (a) The energy spectrum of an 
equilibrium superfluid with T ex —> oo. (b) the indirect exciton 
superfluid with finite, but large r ex (c) The direct exciton with 
short T ex . 



L z — > oo is the normalization length along the z di- 
rection. Obviously, only the photon polarization in the 
plane spanned by k and D k contributes. Note that the 
transition dipole moment D k from the conduction band 
to the valence band at a momentum k is completely dif- 
ferent from the static dipole moment in the dipole-dipole 
interaction Vd(q) in Eqn. [TJ 

In the dilute limit, Vd is relatively weak, so we can 
apply standard Bogoloubov approximation to this sys- 
tem. So in the ESF phase, one can decompose the 
exciton operator into the condensation part and the 
quantum fluctuation part above the condensation bj: — 
\/NS^ + &£. Upto the quadratic terms, the exciton 
Hamiltonian H s j can be diagonalized by Bogoliubov 
transformation: H sf = E(0) +Y,% E $)ptp % where E(0) 

is the condensation energy, E{k) = 



k^k 



2nV d {k)], 

J k ~ U k^k + v %" k * s ^ e Bogoliubov quasi-particle 
annihilation operators with the two coherence factors 



h = 



2E(k) 2- 



The linear term in b k is elim- 



ex 




inated by setting the chemical potential /i = E t 
nVd(0) = {Eg — Eb) + nVd(O). In a stationary state, 
the u. is kept fixed at this value fl3/. As k — ► 0, 
shown in Fig. la where the velocity of the 

Even at 



E(k) 



2ire 2 dn 



quasi-particle is u = y/nV d {0)/M — y M 
T = 0, the number of excitons out of the condensate is: 
n'(T = 0) = | J2k(PJ>z) = J i0?v\ which is the quan- 
turn depletion of the condensate due to the dipole-dipole 
interaction. 

We can decompose the interaction Hamiltonian Hi nt 
in EqnQ]into the coupling to the condensate part H? nt = 



J2k [*ff(^z)(v^ + bo)a kz +h.c] and to the quasi-particle 
part Hf nt — X)fe[*5(^) a fc^ + h.c.]. In the following, we 
study k = and k 7^ respectively. 

3. The photon characteristics at the k = 
mode. The Heisenberg equation of motion of the photon 
annihilation operator is: 



id t (a kz ) = (uj kz - n - i-) (a k J - ig*(k z )VN (2) 
where the average is over the initial zero photon state 



\in) = \BEC)\Q) p h and k is the decay rate of the photon 
due to its coupling to a reservoir to be determined self- 
consistently. The stationary solution of EqnJ5] is [l5| : 



ig*(k z )VN 
(ui kz - n - in/2)' 



(3) 



which is the photon condensation induced by the exciton 
condensation at k = . The photon number distribu- 

where we have 



a\ a kz 



N\g(n/c)\ 2 



tion is n Uk = 

set g*(k z ) around uj kz = fi. The total number of photons 
is N ph = Y,k z n ^kz = N (\g\ 2 D)/ K where D = L z /v g is 
the photon density of states at k = 0. Note that the 
exciton decay rate 70 = \g\ 2 D at k = is independent 
of L z , so is an experimentally measurable quantity. The 
Nph has to be proportional to L z in order to get a finite 
photon density in a given volume L 2 x L z in a stationary 
state. This self-consistency condition sets n = v g /L z — > 
so that 



Nj S(uj kz - fx), N ph = Nj L z /v g 



(4) 



which is independent of L z as required. We showed that 
the power spectrum emitted from the exciton condensate 
has zero width. This conclusion is robust and is inde- 
pendent of any macroscopic details such as how photons 
are coupled to reservoirs. We conclude that the coher- 
ent light emitted from the condensate has the following 
remarkable properties: (1) highly directional: along the 
normal direction (1) highly monochromatic: pinned at 
a single energy given by the chemical potential fi (3) 
high power: proportional to the total number of exci- 
tons. These remarkable properties could be useful to 
build highly powerful opto-electronic device. 

4. Input-Output formalism for a stationary 
state. Now we discuss the properties of emitted pho- 
tons with non-zero in-plane momentum k ^ 0. It is easy 
to see that due to the in-plane momentum conservation, 
the exciton with a fixed in-plane momentum k coupled 
to 3 dimensional photons with the same k, but with dif- 
ferent momenta k z along the z-direction, so we can view 
these photon acting as the bath of the exciton by defin- 
ing = ^ fc g(k)a k . When using the standard input- 
output formalism which treats g( k) non-perturbatively 



under Markov approximation 17[, we find it is conve- 



nient to define the input and output fields as: 



' r(t) = -?^fcr"' (t ' )e 



-i(u k -ix)(t—t ) 



-i(wA,-/i)(t— ii) 



. (5) 



where to — * — 00 and t\ — > 00 are the initial and fi- 
nal time respectively, so to < t < t±. The density of 
states of the photon with a given in-plane momentum 
k is Dr(u!k) = , u>kL . = . It can be shown that the 
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input and output fields obey the Bose commutation rela- 
tions [af (f), a £t(f)] = [a°^{t),a°^{t')]=5 u ,8{t-t>). 

By solving the Heisenberg equation of motion of Eqnfl] 
we find that the output field a|"'(o;) is related to the 
input field by: 

af\u) -- 



[-l + 7]E G n (fc,w + i^)]a^M 



(6) 



normal Green function G n (k,uj) = 
and the anomalous Green function 



where the 

■ u+eg+nVd(k) 
1 u> 2 ~E 2 (k) 

G a (k,u>) — " lVd { k X, . The exciton decay rate in the two 

x ' uj 2 — E 2 (k) ° 

i 1 2 

Green functions are 7^ = D^(fi) \g^(uk = /Lt) which is 
independent of L Z1 so is an experimentally measurable 
quantity. Just from the rotational invariance, we can 
conclude that 7^ ~ const. + \k\ 2 as k — > as shown in 
Fig. lb and lc. Note that the Fourier transformation 
of the Eq. ([5|) leads to uj = u>k — (!>■ Eqn[6] can be 
considered as a S matrix relating the input photon field 
at to — > —00 to the output photon field at t\ — > 00. 

5. Photon number spectrum at k 7^ 0. The an- 
gle resolved power spectrum (ARPS) of the output field 



is S ± $,w) = ^dre-^{a°^{t + r)a^{t) 
inserting EqnJBl one obtains S±(k,u>) = Si(k, uj 



out] , 



lln 2 Vj{k) 

n 2 (Lo) + 1 iE 2 {k) 

k 



By 



(7) 



E 2 (k) + 7g/4. which is shown in 



where il(uj) = ui 
Fig. 2a with different E(k) and 7 e. 

In the strong coupling case k < k*,E(k) < 
7g/2, Si(k,uj) reaches the maximum 7 2 ,n 2 V* ( 2 (fc)/[7 2 ./4 + 



■> 0, E(k) = u\k\ -> 0, 
so the curve has a half 



E 2 (k)] 2 at oj k = fx. As k 
then -> ( J 2+7 . /4)2 , 

width ~ S70 ~ 10~ 4 meV. This is expected, because 
the quasiparticles are not well defined with the decay 
rate 70 much larger than its energy E{k). In the weak 
coupling case E(k) > 7^/2, at the two resonance fre- 
quencies Wfc = /j, ± [E 2 (k) — 7I/4] 1 / 2 , Si(k,w) reaches 
the maximum n 2 V^(k)/E 2 (k). It can be shown that 
when E(k) 3> 7s/2, the width of the two peaks at 
the two resonance frequencies is ~ js, this is expected, 
because the quasi-particle is well defined with energy 
fj, ± [E 2 {k) - 7H/4] 1 / 2 and the half-width 7j j. 

The Momentum Distribution Curve (MDC) S x (k) = 
J2 k , Si(k,w) = j duj k D^(uk)Si(k,uj) is: 



Si(k) 



D^)n 2 V 2 (khj: 
2 [E*(k) + M 



(8) 






FIG. 2: (a) The angle resolved power spectrum (ARPS) 
of the emitted photon with in-plane momentum k. When 
E(k) < ft 7 g/2, there is only one peak with the width 7e/2. 
When E(k) > there are two peaks at the two resonance 

photon frequencies ujk = M± [E 2 (k) — 7 |/4] 1 '' 2 also with width 
7g/2. The Si(k,u) at E(k) = 0.5/ieV and E{k) = 2fj,eV 
are multiplied by 50 and 500 in order to be seen in the fig- 
ure, (b) The zero temperature MDC has a bi-model struc- 
ture consisting of a sharp peak S(k = 0)/L z — Nyo/v g due 
to the condensate at k — superposing on a Lorentzian 
peak with a half width k* ~ W 2 cm~ 1 due to quasi-particle 
excitations at k ^ 0. (c) The zero temperature EDC has 
a bi-model structure consisting of a sharp S function peak 
S(u)/N = yoS(u! — fi) due to the condensate at k = su- 
perposing on a Lorentzian peak with the half envelop width 
7^/2 ~ O.lfieV due to quasi-particle excitations at k ^ 0. 



which is a Lorentian with the half width at k* ~ 
10 _2 cm _1 shown in Fig. 2b. From EqnfJ] one can see 
the condensate contribution at k = is N p h/L z = 
N^o/vg on N, while the contribution from the quasi- 
particle is Si(k -> 0)/L z = 2n l V ^ 0) (x n 2 / 7o . We can 
also calculate the the one photon spacial correlation func- 
tion: 



Gi(r) 



d 2 k 



ik-r 



-k* 1 



(2tt) 2 k 2 + k* 



(9) 



where we can identify the coherence length £ ~ 1/fc* ~ 
40/xra. This coherence length has been measured in Q 
and will be discussed in section 7. 

The Energy Distribution Curve (EDC) Si(w) = 
Y,kSi(k,uj) is: 



Si(w) = Nn x 



d 2 /fc 



7p d 2 (fc) 



(2tt) 2 n 2 (oj)+j 2 E 2 (k) 



NnV${k = 0) M 



47TM 2 



Tfc=o 



(10) 



where f(x) 



1 rw 

x 1-2 



arctg 



} where — tt/2 < 



arctgy < tt/2. From Eqn[H one can see the condensate 
contribution at fc = is n Uk ^/N — r yod(uJk z — fj), while 
the contribution from the quasi-particle Si(lo = 0)/N = 



nVd [ k 'i °' > ■ Both are shown in the Fig. 2c. 



6. Quasi-particle spectrum of a non-equilibrium 
stationary exciton superfluid and its experimental 
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observation. It is important to compare the excitation 
spectra in Fig. 1 . Fig. la is the well know quasi-particle ex- 
citations in an equilibrium superfluid. They are well de- 
fined quasi-particles with infinite lifetime. However, the 
quasi-particles in Fig.lc are not well defined in any length 
scales, because the decay rate is always much larger than 
the energy. Fig. lb is between the two extreme cases. 
When ft < ft*, the quasi-particle is not well defined, the 
ARPS is centered around Uk = ^ with the width 7^. 
The MDC has large values at ft < k*. The EDC has 
large values at u> < 7^ . When k > ft* , the quasi-particles 
is well defined, the ARPS has two well defined quasi- 
particles peaks at wt = /i ± [E 2 (k) — 7|/4] with the 
width 7^. The MDC has very small values at k > k* . 
The EDC has very small values at ui > 7/.. So in the 
long wavelength r > £ ~ 1/ft* ( or small momentum 
k < k* ) limit and long time r > r ex ~ 1/7& (or low 
energy limit oj < 7^ ) limit, there is not a well defined 
superfluid which is consistent with the results achived in 
[1 81 ] . However, in the distance r < £ ~ 1/ft* ( or mo- 
mentum k > k* ) limit and the time r < r ex ~ l/'jk 
(or energy scale u> > 7& ), there is still well defined su- 
perfluid and associated quasi-particle excitations. This is 
the main difference and analogy between the excitation 
spectrum in the equilibrium superfluid in Fig. la and that 
in the non-equilibrium steady state superfluid in Fig. lb. 
Although we derived the Fig. lb from the specific Hamil- 
tonian Eqn.l, we expect it is universal for any station- 
ary pumping-decay system such as exciton-polariton sys- 
tems. Very recently, the elementary excitation spectrum 
of exciton-polariton inside a micro-cavity was measured 
[lfj| and was found to be very similar to that in a helium 4 
superfluid shown in Fig. la except in a small regime near 
ft = 0. We believe this observation is precisely due to 
the excitation spectrum in a non-equilibrium stationary 
superfluid shown in Fig. lb. 

7. Comments on current experiment data and 
possible future experiments. In [J, l5j, the spa- 
tially and spectrally resolved photoluminescence inten- 
sity has a sharp peak at the emitted photon energy 
E = 1.545eV with a width ~ lmeV at the lowest tem- 
perature ~ IK ~ O.lmeV. The energy conservation at 
k z = gives the maximum in-plane momentum k max ~ 
1.545eV/t> fl ~ 4.3 x 10 4 cm _1 where we used the speed of 
the light v g ~ 8.7 x 10 9 cm/s in GaAs 0]. Then the max- 
imum exciton energy E max = uftm W ~ 0.15meV where 
we used the spin wave velocity [J, [5J »~5x 10 5 cm/s. 
The average lifetime of the indirect excitons in the EHBL 
0, HJ is T ex ~ 40ns, then we can estimate the exciton 
decay rate 7fe ~ K/40ns ~ 10~ 4 meV — O.lfieV. At 
the boundary of the two regimes in the Fig. lb where 
E(k*) = uk* = 7fc«/2 = O.lmeV, we can extract ft* = 
2.4 x 10 2 cm -1 <C k max . The typical exciton cloud size 
L ~ 30^tm jij. The number of excitons is N = nL 2 ~ 10 5 
which is comparable to the number of cold atoms inside 
a trap in most cold atom experiments. The central peak 



due to the condensate in the MDC Fig. 2b is broadened 
to fco ~ l/L ~ 10 3 cm _1 which is already larger than the 
half width due to the quasi-particle ft* ~ 10 2 cm -1 . So it 
is impossible to distinguish the bi-model structure in the 
MDC in Fig. 2b at such a small exciton size. The coher- 
ence length was measured in [9]. FromEqnd we find the 
coherence £ ~ 1/ft* ~ 40/im which is slightly larger than 
the exciton cloud size L ~ 30/im. It is easy to see that 
the central peak due to the condensate in the EDC in the 
Fig. 2a is broadened simply due to the change of the lo- 
cal chemical potential from the center to the edge of the 
trap Aw = ^uqL 2 ~ O.lmeV where uo ~ 10~ 12 eVnm~ 2 
0, IH • This value is much larger than the half width due 
to the quasi-particle 7^ ~ 0.1/xeV, so it is impossible to 
distinguish the bi-model structure in the EDC in Fig. 2c 
either. In order to understand if the observed peak is 
indeed due to the exciton condensate, one has to study 
how the bi-model structures shown in Fig. 2 will change 
inside a harmonic trap at a finite temperature ~ IK and 
the effects of both dark and bright excitons. This will be 
discussed in a separate publication. The fine structures 
in the ARPS in the Fig. 2a and phase sensitive homodyne 
measurement 11[ are also needed to test any existence 



of the exciton condensate. 

8. Conclusions We study the power spectrum of 
photons emitted from the exciton superfluid phase in 
semiconductor electron-hole bilayer systems. We find 
that the photons emitted along the direction perpendic- 
ular to the layer are in a coherent state which possesses 
several remarkable properties, while those along all tilted 
directions due to the quasi-particles above the condensate 
show very interesting structures. We determined the an- 
gle resolved power spectrum ( ARPS ), the line shapes 
of both the MDC and the EDC. We also pointed out 
the analogy and difference between the quasi-particles in 
an equilibrium superfluid and those in a non-equilibrium 
stationary superfluid. This difference precisely explained 
the recent experimental observation of excitation spec- 
trum of exciton polariton in a planar microcavity. We 
commented on available experimental data both MDC 
and EDC and also suggested possible future ARPS ex- 
periment to test our theoretical predictions 
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